Non-equilibrium conductance of a three-terminal quantum dot in the Kondo regime: 

Perturbative Renormalization Group 



o : 
o ■ 

<N ■ 

m : 
t> : 

(N . 



c3 



i 

-a 
c 

o 
o 



> 
o 

oo 

On 
O 

in 
o 

-I— > 
c3 



i 

c 

o 
o 



X 



N. Shah and A. Rosch 

Institut fur Theoretische Physik, Universitat zu Koln, 50937 Koln, Germany. 

Motivated by recent experiments, we consider a single-electron transistor in the Kondo regime 
which is coupled to three leads in the presence of large bias voltages. Such a steady-state non- 
equilibrium system is to a large extent governed by a decoherence rate T induced by the current 
through the dot. As the two-terminal conductance turns out to be rather insensitive to F, we study 
the conductance in a three-terminal device using perturbative renormalization group and calculate 
the characteristic splitting of the Kondo resonance. The interplay between potential biases and 
anisotropy in coupling to the three leads determines T and the conditions for strong coupling. 



With the advent of highly tunable strongly correlated 
systems it has not only become possible to systemati- 
cally realize and investigate well-studied theoretical mod- 
els but also to pose new theoretical questions of funda- 
mental nature. Though non-equilibrium physics in itself 
has a long history, a lot remains to be unfolded in the con- 
text of strongly correlated quantum systems. Since the 
single-impurity Kondo model, describing a spin coupled 
to conduction electrons, is not only one of the best un- 
derstood theoretical prototype but also the one that has 
state-of-the-art experimental realization in tunable nano- 
systems involving quantum dots Q, Q , single molecules 
0, 0] and perhaps cold atoms in the near future || , it 
is most natural and crucial to ask: what happens at the 
juncture of Kondo and non-equilibrium physics? 

In equilibrium, the Kondo coupling grows from "weak" 
to "strong" on lowering the temperature and finally the 
spin is screened. An applied bias voltage V influences 
the Kondo physics in at least two ways. First, it pro- 
hibits resonant spin-scattering processes among the leads 
by splitting the relevant Fermi energies by V. More in- 
teresting and less obvious [1,0 HUE] is the fate and 
role of intra-lead resonant spin- flip couplings. Will they 
grow on reducing the energy or will the problem remain 
weakly coupled for large voltages? At least for symmet- 
ric coupling, the Kondo model remains at weak coupling 
for V S> Tk {Tk being the equilibrium Kondo tempera- 
ture), as the large current through the dot induces noise 
and decoheres the resonant spin-flip processes as shown 
m Refs. jEHEEl The physics is then to a large ex- 
tent controlled by the corresponding decoherence rate T. 
Therefore the question arises which experimental quan- 
tity can be used to investigate this second mechanism. 
Unfortunately, the usual two-terminal conductance is al- 
most completely blind to this kind of physics as it it dom- 
inated by (non-resonant) inter- lead m~ocesses|l2j. How- 
ever, a configuration where three leads are coupled to 
a quantum dot has much more flexibility and is indeed 
sensitive to T as shown below. 

Such systems have recently been realized by three 
groups^S Q) E|- Francheschi et al. 0] used a clever 
indirect realization in a dot connected to one lead and 



one mesoscopic wire driven out of equilibrium. The de- 
vice of Leturcq et al., in contrast, is based on a small ring 
with three leads attached jlj|. They managed to operate 
this device in the Kondo regime, where effectively a sin- 
gle spin is localized (or rather delocalized) on the ring. A 
magnetic field allows to fine-tune the wave function and 
therefore to adjust the exchange couplings to the leads. 
A weakly coupled third lead was proposed [lfil Ihil Hsl llflj 
also as a way to experimentally see the splitting of the 
Kondo resonance by a bias voltage, which lies beyond the 
scope of a two-terminal measurement this splitting has 
indeed been observed by these experiments. 

In this letter, we carry out the non-equilibrium per- 
turbative Renormalization Group (RG) analysis ^lj of a 
three-terminal Kondoesque device and calculate the con- 
ductance in the presence of applied voltage biases. We 
model the system by the Hamiltonian 
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where a — 1,2,3 denote the three leads, each with a dis- 
tribution function f a (u>) — 1/(1 + exp/3(cj — /x Q )) with 
the electrochemical potential p a incorporating the ap- 
plied voltages. The key point is that in non-equilibrium, 
all electrons in the energy windows corresponding to the 
applied voltages, contribute to the low energy proper- 
ties and the effective renormalized exchange coupling of 
the electrons from the leads with the local "spin" de- 
pends on their energy. We derive the RG equations 
for these frequency-dependent coupling functions (de- 
noted by g aa ' (cj)) to leading order in the small parameter 
1/ ln[max(|/i Q — fip\ , \u> — fj, a \)/Tic]. The equations read 



dg aa >(uj) 
dhxD 
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with the initial conditions at the bare cutoff D = Dq 
given by g aa ' (lu) — pJ aa > , where p is the density of states. 
The meaning is transparent: At each value of the cut-off 
D, an electron of lead f3 at fip contributes to the renor- 
malization of gu>{D) via a second order virtual process 



2 



only for frequencies that lie within the range of the cut- 
off from pp. By choosing u> = fip, the RG equations are 
first solved for all the constants g a p(np) on the RHS to 
obtain their entire D— dependence, and then in the sec- 
ond step g a a'(fjj) is obtained by using in the RHS the 
pre-evaluated constants at every D. 

As we have emphasized, the crucial element in the non- 
equilibrium is the so-called decoherence rate0 
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cLjg 2 aa ,(Gj)f a (cj)(l-f a ,(Gj)) (3) 



responsible for stopping the RG flow before the strong 
coupling regime given by the above RG equation is 
reached. Implementation of this effect in the RG pro- 
cedure amounts to replacing the O function in Eq. (J2J 

by <d(D— \u> — fip\) (with the short-hand u> = ^/cu 2 + T 2 ) 
and evaluating the value of T self-consistently. 

When the Hamiltonian is assumed to be derived from 
an underlying Anderson model then the J aa > are sepa- 
rable (J QQ ' oc t a t a i where t a is the hybridization) and 
only a single channel couples to the dot with Tk = 
D exp[(3/(2p(Jn+ J22+ J33))] (to one-loop order). Since 
the matrix structure of the general RG equations respects 
this property, we can express g aa > (lu) = 3g u n a n a ' where 
n 2 a = J aa /(Jn + J22 + ^33) is the the lead anisotropy 
parameter (so that n\ + n\ + n 2 = 1) and the initial 
condition is given by g LU (D n ) = 1/[21u(D /Tk)]- 

For concreteness, let us assume that /i 3 i > ^32 > ^21 
with fi a p = \/j, a — fj,p\. When solving the RG equa- 
tions for g a p{pp) as a function of the running cutoff D, 
four different regimes have to be considered. First, for 
D > /i 31 one recovers the usual equilibrium RG equations 
which are easily solved analytically. This is not possi- 
ble in the second regime, /_t 3 i > D > /i 3 2. However, as 
ln[/j 3 i//j 32 ] < In 2 < liifi 3 i/T K (as /j 32 > M21 by defini- 
tion) , the contributions from this intermediate regime are 
small and formally subleading in 1/ ln[/j 3 i/TR-]. We can 
therefore safely take our formulas for D > /j 3 i and use 
them even for D > /j 3 2. In the two remaining regimes, 
/J32 > D > /J21 and finally /i2i > D > again an analyt- 
ical solution can be obtained. When calculating g a p{tu) 
in a second step by simple integration of the RHS of Eq. 
@, we do not use any further approximation and treat 
all four regimes equally (thereby improving also the re- 
sults for gapi/J'p) for /j 3 i > D > /i 3 2 as discussed below). 
The approximation improves with decreasing /J21 and for 
/J21 =0 one recovers correctly the two- lead case for which 
the one loop RG equations can be solved analytically. 

The virtue of having an analytical solution is that it 
provides insight into how the four independent dimen- 
sionless parameters (two Hup/Tk and two rtr a ) enter the 
renormalized g u at various stages of the scaling procedure 
and for what values the strong coupling regime is ap- 
proached. The role of energy scales is exemplified by the 
sheer number of regions (already seven for [i\ < to < /i 3 ) 



into which the entire frequency range gets divided, most 
having a different closed form solution for g^. For the 
frequency range /_t 3 > u> > ^2 + M21, we obtain for D — > 
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while for /i 2 > u> > fJ-i we have 
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As a check for our approximative treatment of the RG 
equations, we compare e.g. 17^3 evaluated by taking 
uj = /j 3 in Eq. with its initial evaluation described 
above (not shown). The difference n\{l/h\ [/J^/Tk] — 
1/ In \p3i/Tk])/2 is even for the worst case (/Z32 = /J31/2) 
only of order In 2/ In 2 [/j 3 i/Tr-] as expected. As Fig. Q] 
shows, there is only a very slight overestimation of g u 
as compared to the straightforward numerical solution of 
Eq. J2J) in this case. In all other cases the numerical and 
analytical curves coincide. 

It is important to check under what conditions the 
problem enters the strong coupling regime where our 
approach is not valid any more. We can investigate 
this question analytically by determining the scale T™ at 
which g fi a diverges in the absence of T. Analyzing Eqs. 
0,JU , for arbitrary \i a p > T* we obtain 



r. = max p; Q ] 

a=l,2,3 



with 



H0a 



(6) 



while if two of the three leads have the same electrochem- 
ical potential (more precisely /ii2 < T», say) then 



T* = 
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FIG. 1: Upper panel shows the scaling of as the cut-off is 
decreased down to zero from D — 300Tk, for fii — IOOTjc, 
fi2 = 150Tk, = 350Tk and ni = 112 = n^. The light 
grey plot shows the g^ obtained if Y were to be T*. Lower 
panel shows the renormalized gu{D — 0) for /ii = IOOTr-, 
jj,2 = 225Tk, fJ.3 ~ 350Tx for different ni,n2,nz shown in 
the legends. In the same line style as numerical plots, also 
shown in both panels are the plots for g u obtained analytically 
(upper curves when they do not coincide). 



For r = T» at least one of the coupling constants diverges 
oc l/ln[r/T*] (this is how we determine T"; note that 
x = \/ x 2 + r 2 in the terms ln[x/Tff]). The perturbative 
RG is therefore controlled by the ratio T/T*. The larger 
is the ratio as compared to one, the further away we are 
from the strong coupling coupling regime, and the more 
accurate is our analysis. 

The upper panel of Fig. ^ is a visual display of the 
scaling process: g u for a generic case is plotted as the 
cut-off D is gradually decreased; one can see the role of 
energy scales /X21, /i3i and D, all measured in units of 
Tk- To clarify the role of T, we also plot (in light grey) 
the g w obtained if T were to be T*. In the lower panel 
we illustrate how the interplay between lead anisotropics 
and potential differences determines the renormalization 
of g u . For example, when one lead a is very weakly 
coupled, there is no feature at u> = fi a . 

The physical quantity of interest is the non-equilibrium 
conductance matrix, which, within the precision of our 
approach, can be evaluated by plugging the renormal- 
ized coupling constants in the second order perturbative 
formula to obtain 
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(8) 



where we have restricted to zero temperature and intro- 
duced the quantum of conductance, G° = e 2 /irh. Let 
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FIG. 2: Conductance Gn in the "two-lead spectroscopy" case 
(n\ — 0.002) for ^32 = IOOTk and different values of n| and 
n§. Inset shows F and T* on a log scale, as a function of n\. 



us first consider the case when lead one is a "prob- 
ing lead" i.e. it is very weakly coupled to the dot 
(n\ <C 1). Then Gn(/Lti) is simply given by the first 
term of Eq. (JBJ (the second term being negligible ow- 
ing to very weak dependence of g 2 {e) on /ii) and traces 
the two-lead spectral function of the quantum dot; we 
have explicitly confirmed this equivalence. If we as- 
sume an underlying Anderson model as above, we obtain 
G n (fii)/G° = 3tt 2 7i 2 (1 - nDg 2 ^ with g^ given by Eq. 
(|5|). As /ii is swept, Gn peaks up at ^2 and at 113 with 
the peak height at ^2 given by 
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while that at 113 obtained by interchanging n 2 <-> n 3 in 
the above expression. For n\ <C 1 the formula for T* 
(Eq. ©) reduces to Eq. Q which now holds not only 
for /ii = H2 but for all values of n\. 

In Fig. El we plot the conductance for different values 
of the relative anisotropy between couplings to leads two 
and three, for a given value of potential difference /Z32, 
to emphasize the role of anisotropy. The current from 
lead two to three is largest for n 2 ~ n 3 = 0.5 and accord- 
ingly r (same for all /zi as n\ <C 1) is maximum at this 
point . On the other hand, T* is minimum at this point 
and the ratio of T/T* largest (see inset). In the extreme 
anisotropy case when only a single lead a is strongly cou- 
pled to the dot, one will measure the equilibrium spectral 
function (g(oj) is strongly renormalized only at that fi a ). 
In Fig. |3 we show how the conductance changes as the 
potential difference /132 is increased, for a given relative 
anisotropy. To clarify how T determines the rounding of 
the peaks, the light line shows Gn if T were to be T* (for 
^32 = lOOTx). With increasing /x 32 , the value of T gets 
larger (since the current increases ) and T* gets smaller, 
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FIG. 3: Conductance Gii in the "two-lead spectroscopy" case 
(n? = 0.002) for n 2 = 0.599, n 3 = 0.399 and different values 
of (j,32- The light curve shows Gii if F were to be T*. The 
inset shows T and T* on a log scale, as a function /i32. 



making the ratio r/T* increasingly larger (see inset). 

To study other cases, we phrase the question: how does 
the conductance change as we increase the coupling to 
lead one? We hold the relative anisotropy between leads 
two and three to be constant as we increase the value of 
7ii. As can be seen from Fig. as we start increasing 
ni, the overall conductance Gu begins to grow. Further, 
the second term in Eq. JSJl is not negligible as in the 
"two-lead spectroscopy" case and makes a negative con- 
tribution; the inset plots the peak height Gn{m = (13) 
as well as the contribution it gets from the first term to 
make the increasing deviation between the two evident. 
As we keep increasing m, we start approaching another 
limit in which lead one is strongly coupled and leads two 
and three weakly coupled to the quantum dot: the peaks 
continue to grow, however, the overall conductance for 
other values of /ii begins to drop down. For the sake of 
presentation we have chosen to focus on Gn although it 
should be clear that analogous arguments hold for the 
remaining elements of the conductance matrix. 

As can be seen from the insets of Figs. I2I3I4I the ratio 
r/T* can become quite large for some parameter range 
making our perturbativc RG treatment extremely accu- 
rate. Further more, only when the system approaches 
equilibrium and T* — > Tk, do we find that r — > T* 
and the strong coupling regime is approached. This 
happens either when /X32 < Tk for the n\ <C 1 case 
and all fj, a p < Tk for the general case, or in the ex- 
treme anisotropy limit when only one of the three leads 
is strongly coupled to the dot irrespective of the values of 
fi a p. Our finding is in accordance with the expectation 
for a Kondo model derived from an Anderson model: a 
strong coupling fixed point exists only in equilibrium [3]. 

Possible extensions to this work include studying the 
case of a generalized Kondo model with no underlying 
Anderson model and considering the effect of excited lev- 
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FIG. 4: Plots of Gn for increasing values of n\ (ordered by 
increasing peak heights) for fixed = 1.5 and /132 = 

IOOTk . The corner inset shows Y and T» on a log scale, as 
a function of nf for fit — fi2, fJ-3- Other inset shows the peak 
height Gn(/ii = H3) as well as the part coming only from the 
first term of Eq. [H] (dashed line) plotted against n\ . 



els or of a Zeeman field (B 3> Tk)- For the specific 
set-u p o f a quantum ring geometry enclosing a magnetic 
flux hj it might be interesting to investigate the chiral 
component to the current [20j in the Kondo regime and 
in non-equilibrium. 
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